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We have investigated the effect of imperfect circular polarization on the angle-resolved photoe¬ 
mission spectroscopy signal, using graphene as a prototypical system that can be understood within 
tight-binding formalism. We found that perfect left- and right-circularly polarized lights give the 
same photoelectron intensity distribution around a constant energy contour of the graphene tt band. 

On the other hand, upon breaking the purity of the polarization, photoelectron intensity starts to 
show circular dichroism, which is enhanced with further increasing the imperfection. Our results 
predict the existence of an additional factor for the circular dichroism observed in the photoemission 
signal from graphene and hence suggest the importance of experimental conditions to understand 
circular dichroism observed via photoemission spectroscopy. 


I. INTRODUCTION 

Circular dichroism has been one of the powerful 
methodologies to extract information on spin and/or 
orbital properties of charge carriers in solid state sys¬ 
tems [iMsI. In addition, recent study on graphene has 
suggested that information on Berry phase can also be 
obtained via the circular dichroism Q , extending a previ¬ 
ous approach of the direct measurement of Berry phase 
using linearly polarized lights for the same system [5|. 
These results provide an experimental evidence that the 
quantum mechanical phases can be probed by photoemis¬ 
sion spectroscopy, previously not believed to be possible, 
and hence constitute the first band specific measurements 
of Berry phase. 

These interesting observations have been possible due 
to the simple geometric structure of graphene, allowing 
us to obtain the explicit form of the initial electronic 
states within the tight-binding formalism [6|. Shirley et 
al, 0 calculated photoelectron intensity, which is the 
absolute square of the transition matrix element Mk = 
(/k| l^/k), where j^/k) is a tight-binding eigen¬ 

state, l/k) is a plane-wave final state, and 77^^^ = A • p. 
Here, A is a light polarization and p = —ihV is the mo¬ 
mentum operator, where h is the Planck’s constant. 

This approach reproduces the photoelectron intensity 
for the linearly polarized light along the x-axis denoted 
in Fig. 1(a), i. e., X-polarization (A = A^^x), whereas it 
is not successful in reproducing full polarization depen¬ 
dence of the photoelectron intensity, e. g., when the light 
polarization is rotated by 90 degrees, i. e., Y-polarization 
(A = Ayy) 0. This issue originates from the applica¬ 
tion of p = —i W to the tight-binding eigenstates. When 
the tight-binding Hamiltonian is intrinsically non-local, 
the derivative in real space for the tight-binding eigen¬ 
state does not work [H, ^ . This suggests that the agree¬ 
ment even for X-polarization using p Q could be fortu¬ 
itous. In order to solve this issue, Hwang et al, [5[ intro¬ 
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duced an alternative approach replacing the derivative 
to the commutation relation, i. e., p/mo = —iW/mo, 
where mo is free-electron mass, to v = [r, /ih, where 
r = i^(Vk,9/eJ in the k-representation p. [iQj. 

To avoid this issue, Liu et al, have assumed real¬ 
istic final states, i. e., the Block sum of the Wannier 
states, while still using p Q. With this setup, they have 
claimed that the observed circular dichroism originates 
from Berry phase of graphene, which is challenged by 
a study on the circular dichroism as a function of elec¬ 
tron binding energy 0. In fact, another experimen¬ 
tal study shows that the circular dichroism varies upon 
changing photon energy, emphasizing the role of the final 
state effect in understanding photoemission signal from 
graphene [T2|. These controversies give rise to a funda¬ 
mental question on the origin of the circular dichroism. 

Here we report calculated photoelectron intensity of 
graphene for circularly polarized light. The photoemis¬ 
sion matrix element was constructed using the velocity 
operator 0i within the tight-binding formalism. We 
found that the photoelectron intensity for left- and right- 
circularly polarized light (LCP and RCP, respectively) 
does not show any difference for perfect circular polar- 



FIG. 1. (a) Schematic of graphene in real space, (b) The 
first Brillouin zone of graphene. Here, bi = 5(0,1), b 2 = 

5 — 2^ and b3=5^^,—2^ are the three vectors con¬ 

necting the in-plane nearest neighbor atoms where b = 1.42 A. 
The positions of the K and K' points are (|^, O) and , O), 
respectively, where a = \/35. 
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ization. On the other hand, upon decreasing the purity 
of the polarization, the photoelectron intensity starts to 
exhibit circular dichroism. Our results indicate that an 
experimental condition in conjunction with the intrinsic 
effects such as the isospin of charge carriers and the final 
states [ill plays an important role in determining photo¬ 
electron intensity. 


II. TIGHT-BINDING FORMALISM 


In the presence of the vector potential A, the Hamil¬ 
tonian is obtained by Peierls substitution, i. e., k —> k — 
^ A. Then, the interaction Hamiltonian is obtained 
by the first-order term of A from (k — -^A) — i^^(k). 


and represented as — ^A • v using the velocity operator 




/ih, where f = Vk in the k-representation. 


h is the Planck’s constant, e is the charge of an electron, 
c is the speed of light. Then the interaction Hamiltonian 
becomes [5[: 


Figure 1(a) shows the geometric structure of graphene 
with two carbon sublattices, A and B. In order to de¬ 
scribe the electron energy eigenvalues and wavefunctions, 
we have used tight-binding formalism for the Pz orbital 
of each sublattice using the in-plane nearest-neighbor 
(A-B) hopping integral, to (we confine our interest to 
this single parameter for simplicity), which corresponds 
to — 7 q in the well-known Slonczewski-Weiss-McClure 
model [a Because our basis set for the tight- 

binding model has more than one non-equivalent or¬ 
bitals, there exist two different hopping parameters dif¬ 
fering only in the sign that give exactly the same electron 
band structure. In this sense, the absolute magnitude of 
inter-orbital hopping integrals within the empirical tight- 
binding Hamiltonian Htb between non-equivalent states 
(i. e., ( 011 Htb 102 ), where the localized orbitals |0i) and 
102 ) are not equivalent) in any material have only been 
speculated theoretically, whereas it has recently been 
proved that the sign of the hopping integrals for both 
single- and double-layer graphene can be experimentally 
determined uniquely using the angle-resolved photoemis¬ 
sion spectroscopy (ARPES) technique We have used 
001 = 3.16 eV and to < 0, the values in Table H of 
Griineis et al, [lil and in Fig. 8 of Hwang et al. 0, 
respectively. 

The tight-binding Hamiltonian of graphene for two- 
dimensional wavevector k is as follows [^, [6[ : 

/ 0 tog{k) \ 

= , ( 1 ) 
\tog*{k) 0 J 

using a basis set composed of Bloch sums of localized 
orbitals on each sublattice: 

3 

5(k) = exp(ik • hi) (2) 

i=l 


with b0s defined as in Fig. 1(a), and 


1 
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1 
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( 3 ) 


FF“‘(k) = - —A 

he 


0 fo Vk5(k) 

foVk5*(k) 0 


( 5 ) 


Note here that because g{k) depends only on the and 
ky components of the wavevector k, there is no contribu¬ 
tion arising from the z component of the vector poten¬ 
tial Az within this tight-binding model. In real measure¬ 
ments, the light with a nonzero polarization component 
along the 2 ; direction will give rise to an additive isotropic 
term to the photoelectron intensity that is independent 
of the in-plane polarization of the light. 


III. THE ANALYSIS OF THE 
PHOTOELECTRON INTENSITY 

The photoelectron intensity is described by the ab¬ 
solute square of the transition matrix element M^k = 
(/k| i^^^^(k) l^sk), where I'lps'k) is graphene eigenstate 
with the band index s = ±1 for conduction (+) and 
valence (—) bands, and |/k) is the plane-wave final state 
projected onto the Pz orbitals of graphene. For graphene, 
we may use 


Here, we neglect the kz dependence of the final state, 
as done in the previous works 0 0 0. For a few 
tens eV photons, typical for ARPES measurements for 
graphene [3, H, [lil, [ill, El, EH , kz of the plane-wave fi¬ 
nal state is much larger than k^ and ky , leading to only a 
small variation in kz with any change in kx and ky . In ad¬ 
dition, the study on i. e., photon energy dependence, 
is beyond the capability of the tight-binding formalism, 
but can be achieved through the first principles calcula¬ 
tions EH. 

When we consider the case where k is very close to 
the Dirac point K as denoted in Fig. 1(b), and define 
q = k — K (|q| ^ |K|), Eq. ([2]) becomes 

g{q + K)K-^b{q:,-iqy) , (7) 

and the Hamiltonians become 

i7°(q + K) « —^ (g® CTj; + qy ay) , 


0 

1 


k 


( 4 ) 


(8) 
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and 

i/-‘(q + K)«|^6to(A-a) , (9) 

where cr is the Pauli matrix. Then the electron energy 
eigenvalues and wavefunctions of the are given by 
-Bsk = |to| S |q| and 

1 / g-i0q/2 \ 

respectively, when is the angle between q and the -\-kx 
direction. With this setup, the photoemission matrix 
elements for X- and Y-polarizations are given by 

~ exp(-i6'q/2) + sexp(0q/2), (11) 

and 

-^+rk°* ~ exp(-i0q/2) - sexp(6»q/2), (12) 

respectively, of which absolute square, i. e., Igk = 
\Mskf^ represents photoelectron intensity that repro¬ 
duces the previous experimental results for linearly po¬ 
larized lights ii- Resultantly, the variation of photo¬ 
electron intensity depending on the polarity of light re¬ 
veals the pseudospin nature of charge carriers in single- 
and double-layer graphene, i. e., Berry phase of tt and 
27r, respectively, and the signs of hopping integral in the 
tight-binding Hamiltonian for double-layer graphene or 
graphite [5|. 

In the same analogy, for the circularly polarized light, 
i. e., A = AxX ^ the photoemission matrix ele¬ 

ment is given by 

Afsk = ^ (5 exp(i6>q/2) + exp(-i6>q/2)) 

± ^{s exp(i6>q/2) - exp(-i6>q/2)) , (13) 

where + and — correspond to LCP and RCP, respec¬ 
tively. For the states above the Dirac energy, i. e., s = -hi 

M+ik = cos(6>q/2) + i Ay sin(0q/2), (14) 

and 

M+il = A:c cos(6>q/2) - i Ay sin(0q/2), (15) 

and for the states below the Dirac energy, i. e., 5 = — 1 

Aa, sin(6>q/2) + cos(6>q/2), (16) 

and 

= —i Aa, sin(6>q/2) — A^ cos(6>q/2). (17) 

It follows that 

J+ik = Al cos^(6'q/2) + Al sin^(6'q/2), (18) 


and 

/_i k = A^ sin^ (^q/2) + Ay cos^ (6>q/2), (19) 

regardless of the chirality of the light. This simple al¬ 
gebra indicates that, within the tight-binding formalism 
for |q| ^ |K|, the photoelectron intensity of graphene 
does not show circular dichroism. Moreover, for per¬ 
fect circular polarization, for which \Ax\ = |A^|, 
is isotropic around the constant energy contour. How¬ 
ever, in real measurements, the photoelectron intensity 
for X- and Y-polarization geometries are different due 
to nonzero polarization component along the z direction, 
e. g., Ax exhibits a finite z component when projected 
to the sample surface, whereas Ay has negligible out-of- 
plane component [5[ , when the photoelectron intensity is 
closely related to the scattering probability of the real Pz 
orbitals for the out-of-plane component of light polariza¬ 
tion, which is beyond the capability of the tight-binding 
approach. For example, at a photon energy of 50 eV, the 
ratio of ^21.4 [5[ that we have used 

throughout our study. This ratio is not controllable, but 
determined for each experimental setup for each photon 
energy. This is also applied to the circular polarization, 
for which A = Ax x:t i Ay y is projected onto the sample 
surface. 

These theoretical results are summarized in Fig. 2 at 
several different energies with respect to the Dirac en¬ 
ergy, Ed. For LCP (Fig. 2(a)), the constant energy con¬ 
tour shows a crescent-like shape with minimum intensity 
in the first Brillouin zone (1st B.Z.) at — Ed = 0.4 eV 
and a point-like constant energy map at Ed- The inten¬ 
sity distribution is reversed at — F^d = —0.4 eV, with 
respect to F^d, with maximum intensity in the 1st B.Z. 
While the overall shape shows the characteristic conical 
dispersion of graphene, the intensiW distribution is sim¬ 
ilar to the case of X-polarization [3, Q . It is important 


(a) (b) 



- 0.15 0.00 0.15 - 0.15 0.00 0.15 


k^-K (A-i) 


FIG. 2. (a, b) Calculated photoelectron intensity of graphene 
for (a) left- and (b) right- circularly polarized lights. An ar¬ 
bitrary energy broadening of 0.10 eV has been used to quali¬ 
tatively compare with experimental results BUI- 
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to note that, for linear polarizations, the crescent-like 
shape is understood by the interference between photo¬ 
electrons emitted from two carbon sublattices Q. On 
the other hand, for circular polarizations, it originates 
from the relatively weak photoelectron intensity for Y- 
polarization as discussed above. As a result, the differ¬ 
ence in photoelectron intensity between LCP and RCP 
is not determined by the properties of initial states, but 
given by the experimental geometry, and hence we ex¬ 
pect the same intensity distribution for RCP as shown in 
Fig. 2(b). 

Note also that our results (Eqs. (14-15) and (18)) for 
the conduction band are different from those in the pre¬ 
vious study (Eq. (3) of the Ref. Q): 


1 = 




( 20 ) 

where -|- and — signs correspond to LCP and RCP, re¬ 
spectively. In this study Q, they have claimed that the 
dipole transition matrix elements and for the x and 
y components of the vector potential have a relation of 
~ at 30 eV [i3- This setup results in following pho¬ 
toelectron intensity (Eq. (4) of Ref. Q) for the perfect 
circular polarization, i. e., \Ax\ = \Ay\ = |A|, 


J = 4 A^ |cos(6»/2) ± sin(6»/2)|^ . (21) 




0.15 


0.00 


-0.15 


(a) 


Polarization purity: 100 % 



Polarization purity: 80 % 


E-Eo = 0A eV 


(c) Polarization purity: 50 % 


Max. 

Min. - 

-0.15 0.00 0.15 -0.15 0.00 0.15 





The difference between the two theoretical approaches 
(Eq. ([18]) vs. Eq. ([2T]) ) originates from the evolution of 
the tight-binding Hamiltonian. Eor example, for |q| ^ 
|K|, the Hamiltonian is expressed via Pauli matrices, 
i. e., H = • q (Eq. (8)), where is the Eermi ve¬ 

locity, with which we can obtain the tight-binding eigen¬ 
states (used in both works done by us (Eq. (10)) and 
Liu et al. 0). Hence, the Peierls substitution natu¬ 
rally leads to = a • A both of which compo¬ 

nents exhibit a strong influence on the photoemission 
matrix element when applied to the spinor eigenstate 
of graphene ^ as we have done in our study. On the 
other hand, Liu et al. Q, have applied the local Hamil¬ 
tonian, H = -p‘^12m + E(r), to the tight-binding eigen¬ 
states (Eq. (10)) obtained by the nonlocal tight-binding 
Hamiltonian (Eq. (1)) [9|. Consequently, the y compo¬ 
nent of the matrix element differs by the imaginary num¬ 
ber “i” arising from Gy (compare Eqs. (14-15) and (|2T]) 1 
resulting in the completely different photoelectron inten¬ 
sity distributions [this issue is eliminated when the realis¬ 
tic eigenstates, e. g., maximally localized Wannier func¬ 
tions obtained by ab initio calculations [lH, are used in 
conjunction with the local Hamiltonian]. 

In real measurements, in contrast to our prediction, 
the photoelectron intensity from graphene exhibits cir¬ 
cular dichroism B ini [l2j . The intensity maximum 
around a constant energy contour rotates by 180 de¬ 
grees upon changing the chirality of light with an en¬ 
ergy of 30 eV 0. Such a dichroic effect varies at dif¬ 
ferent photon energies M , which has been attributed 
to the symmetry of the final states (d-like partial waves 


^x-K(A-4 

FIG. 3. (a-c) Calculated photoelectron intensity of graphene 
for LCP (left panels) and RCP (right panels) at E — Ed = 
0.4 eV for 100 %, 80 %, and 50 % circular polarization. An 
arbitrary energy broadening of 0.10 eV has been used. The 
photon energy is assumed to be 50 eV and the angle of inci¬ 
dent photons and surface normal is 26 degrees, adapting the 
experimental geometry in Ref. @]. 


above a photon energy of 52 eV and s or p-like partial 
waves below 52 eV [l9|). The photon energy dependence 
suggests that, when the final state effect is minimized 
around 52 eV, the still observed circular dichroism orig¬ 
inates from the psuedospin nature of charge carriers in 
graphene [T2|. Another study shows that, at a similar 
photon energy of 50 eV, the circular dichroism changes 
as a function of electron binding energy, which suggests 
a possibility of many-body interactions as an origin of 
the observed dichroism [ll|. These ambiguities in un¬ 
derstanding experimental results may suggest the exis¬ 
tence of an extrinsic factor contributing to the circular 
dichroism that changes at different experimental geome¬ 
tries and photon energies, such as imperfect circular po¬ 
larization. In fact, the polarization purity for LCP and 
RCP in the experimental work for the Berry phase sce¬ 
nario is ^80 % 0. 

In Eig. 3, we show the effect of the imperfect polariza¬ 
tion on the rotation of intensity maximum upon chang- 












5 


ing the chirality of light. We assume that the angle 
between the photon incident and the electron detector 
angles is 55 degrees (e. gj, the experimental setup for 
the previous experiment [5|). Then the angle between 
the photon incident and the sample normal angles is de¬ 
termined by photon energies. For example, at 50 eV, 
where the K points is tilted by ^29 degrees with re¬ 
spect to the sample normal, the angle between the pho¬ 
ton incident and the sample normal angles is ^26 de¬ 
grees (14 and 31 degrees for 30 and 70 eV photons, 
respectively). In addition, for imperfect circular polar¬ 
ization, the major axis of the elliptical polarization is 
tilted by ±45 degrees for LCP and RCP, respectively, 
with respect to -\-x direction. As a result, the elliptical 
polarization Ax x ± iAyy projected to the sample sur¬ 
face becomes Ax cos 26° cos 45° i Ay cos 26° sin 45° and 
—Aa,sin45° ± i Ay cos 45° along x- and ^/-axis, respec¬ 
tively, where (1 — (a — 6)/(a ± b)) x 100 % defines the 
purity of the polarization. When the purity of the po¬ 
larization is 100 % (Fig. 3(a)), the intensity maximum 
around the constant energy contour at E — =0.4 eV 

stays the same upon changing the chirality of light as 
denoted by the white arrows. When the purity of the 
polarization decreases, the intensity maximum becomes 
separated as denoted by the white arrows in Fig. 3(b) 
giving rise to the circular dichroism. The circular dichro- 
ism becomes stronger with further decreasing purity of 
the polarization as shown in Fig. 3(c). 

To further discuss the discrepancy of our results com¬ 
pared to experimental results, we directly compare the 
rotational angles in both our and previous studies. When 
the final state effect is minimal at 52 eV M , i. e., in¬ 
trinsic properties of graphene dominate photoelectron in¬ 
tensities ^_the rotation of intensity maxima is only ^40 
degrees [ 13 , different from the rotation by 180 degrees 
corresponding to the Berry phase scenario Q. While 
the tight-binding approach that we have used excellently 
reproduces the experimental data for linearly polarized 
light with an energy of 50 eV the observed (^40 de¬ 
grees) and predicted (^60 degrees via the first principles 
calculations) dichroism 0 is comparable to our result 
of 40-60 degrees for the polarization purity of 80-50 % 
(Fig. 3). 

Now let us restrict our discussion to the case where the 
photon energy only changes photon polarization with re¬ 
spect to the sample surface. With this set up, the simula¬ 
tions for left-circularly polarized light with a polarization 
vector corresponding several photon energies (30, 50, and 
70 eV in Figs. 4(a-c), respectively). Here, the angle be¬ 
tween incident photons and the analyzer is assumed to 
be 55 degrees. This simulation does not show a notable 
rotation of photoelectron intensity around a constant en¬ 
ergy contour upon changing photon energy, suggesting 
that the final state effect [13 plays a dominant role in 
reproducing the photon energy dependence, while the 


experimental condition in conjunction with intrinsic ef¬ 
fects, such as the pseudopin nature of charge carriers 0 
and the many-body effects 0 , also plays a finite role in 
rotating photoelectron intensity upon changing the chi¬ 
rality of light at constant photon energy. Overall, the 
theoretical approach within tight-binding formalism pre¬ 
dicts an additional factor for the circular dichroism and 
hence suggests the importance of experimental conditions 
to understand circular dichroism observed via photoemis¬ 
sion spectroscopy. 


(a) Purity: 50 % (b) Purity: 50 % (c) Purity: 50 % 
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FIG. 4. (a-c) Calculated photoelectron intensity of graphene 
for RCP at E — Er, = 0.4 eV. An arbitrary energy broadening 
of 0.10 eV has been used. The photon energy is assumed to 
be (a) 30, (b) 50, and (c) 70 eV, where the angle of incident 
photons and surface normal is 14, 26, and 31 degrees, adapting 
the experimental geometry in Ref. [g]- 


IV. CONCLUSION 

We have investigated the role of imperfect circular po¬ 
larization in the circular dichroism of photoelectron in¬ 
tensity observed in graphene. Within the tight-binding 
formalism, we found that the calculated photoemission 
matrix element does not predict any difference between 
left- and right-circular polarization. However, we found 
that circular dichroism is developed and enhanced with 
decreasing purity of the polarization. Our results implies 
that the experimental conditions should be taken into 
account to understand circular dichroism, which invites 
further experimental and theoretical investigation to un¬ 
derstand the origin of the observed circular dichroism in 
graphene. 
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